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ABSTRACT 

Extending previous studies of nonthermal electron transport in solar flares 
which include the effects of collisional energy diffusion and thermalization of fast 
electrons, we present an analytic method to infer more accurate estimates of the 
accelerated electron spectrum in solar flares from observations of the hard X-ray 
spectrum. Unlike for the standard cold-target model, the spatial characteristics 
of the flaring region, especially the necessity to consider a finite volume of hot 
plasma in the source, need to be taken into account in order to correctly obtain 
the injected electron spectrum from the source-integrated electron flux spectrum 
(a quantity straightforwardly obtained from hard X-ray observations). We show 
that the effect of electron thermalization can be significant enough to nullify the 
need to introduce an ad hoc low-energy cutoff to the injected electron spectrum in 
order to keep the injected power in non-thermal electrons at a reasonable value. 
Rather the suppression of the inferred low-energy end of the injected spectrum 
compared to that deduced from a cold-target analysis allows the inference from 
hard X-ray observations of a more realistic energy in injected non-thermal elec¬ 
trons in solar flares. 

Subject headings: Sun: corona - Sun: flares - Sun: X-rays, gamma rays 
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Introduction 


The copious hard X-ray emission produced during solar flares is the main evidence 


Holman et al. 

2011 ; 

Kontar et al. 

2011 ) 


20111) . These accelerated electrons propagate through the 


surrounding plasma, where electron-ion collisions give rise to the observed bremsstrahlung 


hard X -rays . The collisions with ambient particles, principally other electrons (e.g., 


Brown 


1972 


Emslie 


19781 ). are responsible for energy transfer from the accelerated particles to the 


ambient plasma. 


In a cold-target approximation, where the 


are dominated by systematic energy loss (e.g., 


dynamics of hard X-ray emitting electrons 


Brown 


19711) . only a small fraction (~ 10 5 ) 


of the energy in the accelerated electrons is emitted as bremsstrahlung hard X-rays. An 
observed hard X-ray flux thus translates into a much higher accelerated electron energy 
content, to the extent that the energy in such accelerated electrons can be comparable to 


the flare (see 


Emslie et al. 

to 

0 

0 

fefe_ 

2005, 

2012 ) 


Solar flare hard X-ray spectra /(e) (photons cm 2 s 1 keV 1 at the Earth) in the 


nonthermal domain are typically quit e steep, w ith power-law forms /(e 


indices 7 ~ (3 — 6 ) (see, e.g., 


Kontar et ah 


2011 


r~7 


and spectral 


, for a review). Using a cold target model 


(that retains only the effect of energy loss in the dynamics of emitting electrons) requires 
that the injected electron flux spectra F 0 (E 0 ) oc E£ 5 (el ectrons cm -2 s -1 keV^ 1 ) are 


similarly steep, with a power-law indices <5 = 7 + 1 (e.g., 


Brown 


19711). Since the total 


injected energy flux £ = J 0 °° E 0 F 0 (E 0 ) dE 0 diverges at the lower limit for such steep power 
laws, the concept of a “low-energy cutoff” E c is frequently assumed in or der to keep t he 


value of £ = J^° E 0 F 0 (E 0 ) dE 0 finite. Indeed the value of E c is usually flHolman et ah 


20031 ) taken to be the maximum value consistent with the hard X-ray data, resulting in the 
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minimum energy in the nonthermal electrons that is consistent^] with the data (see, e.g., 


Holman et al. 


2003, 


201 ill 


,0 


Observations indicate that fast electrons are accelerated in, and subsequently move 
through, the hot (~10' K) flaring plasma in the corona. In some cases, the coronal density 


is sufficiently high to arrest the electrons wholly within t 


thick-target coronal source (e.g., 


Veronig Brown 


2004; 


ris coronal reg ion, resulting in a 


Xu ct al. 


2008 1. In other cases, 


the coronal density is sufficiently low that the electrons emerge, with somewhat reduced 
energy, from the coronal region and then impact the relatively cool ( 10 4 K) gas in the 
chromosphere, producing the commonly observed “footpoint-dominated” flares. The plasma 
in the corona and in the preflare chromosphere is rapidly heated to temperatures in excess 
of ~10' K, producing the commonly-observed soft X-ray emission over the extent of the 
flaring loop. This heating (or direct heating in magnetic energy release) necessarily renders 
a significant portion of the target “warm”; i.e., such that, for an appreciable portion of the 
injected electron population, the injected energy E 0 is comparable to the thermal energy 
kT, where k is Boltzmann’s constant and T the target temperature. 


Ems lie fl2003f l included consideration of the finite temperature of the target in 


modifying the systematic energy loss rate of the accelerated electrons; such considerations 
come into play as the electron energies E approach a few kT. He showed that this resulted 
in a lowering of the requir ed energy flux £ of injected electrons for a prescribed hard X-ray 


intensity. Further, 


G all oway et al. 


(120051 ) have emphasized the role of energy diffusion, 


which is also important at energies of a few kT and is a necessary in gredie nt for describing 


thermalization of the fast electrons in a warm target. 


Jeffrey et al. 


(20141) showed that 


4 Because of the dominance of the (even steeper) thermal hard X-ray component at low 
photon energies e, lower values of E c are permitted by the data but result in unjustihbaly 
large values of the injected energy flux £ = E 0 F 0 (E 0 ) dE 0 . 
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th e transport o f elect rons is more complicated than assumed by either Em si id ([20031) or 


Gallowa y et al. 


(120051 ). inasmuch as the effects of diffusion in both energy and space must 


be included in a self-consistent analysis of electron transport in a warm target. 


In this paper we follow 


Jeffrey et al . 


Gallowa y et al. 


(120051 ). 


Goncharov et al. 


(120101 ). and 


(20141). We show that thermalization of fast electrons in a warm ambient 


target significantly changes the evolution of the electron energy spectrum comp ared 


in a finite temperature target analysis that neglects diffusion in energy 


even more so compared to the standard cold target model (e.g., 


Brown 


Imslie 


o that 


2003]), and 


19711). We also 


emphasize that the appropriate treatment of the thermalization of electrons injected into a 
warm target results into a higher hard X-ray yield per electron, so that significantly fewer 
supra-thermal accelerated electrons are required to be injected into the target in order to 
produce a given observed hard X-ray flux. Moreover, our analysis shows that, contrary to 
the case of a cold target (in which the spatially-integrated hard X-ray yield is independent 
of the density profile of the target), in a warm target one does need to take into account the 
spatial characteristics of the emitting region, in particular the extent of the warm target 
compared to that of the overall flaring region. 


The outline of the paper is as follows. First, we quantitatively evaluate the effects 
of thermalization on electron transport, using both analytical (Sections 12.11 and 12.21) 
and numerical (Section I2.3[) methods, which we compare in Section 12.41 We find that 
warm-target thermalization of electrons considerably reduces the flux of injected electrons 
compared to that obtained in a cold-target model that assumes a low-energy 
the thermalization limit derived here, and even compared to that in a model 
that includes warm-target effects on the secular energy loss rate but neglects diffusion 
in energy. In Section |3] we summarize the results and point out that the magnitude 
of the effect is sufficiently large so that a low-energy cutoff E c in the injected power 


cutoff 


oelow 


(Emslie 


2003) 
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V = A f E 0 E 0 (E 0 ) dE 0 (where A (cm 2 ) is the injection area) is a natural consequence, 
rather than an ad hoc assumption adopted a posteriori. We therefore urge discontinuance 
of cold target modeling for a warm plasma and its resultant need to impose an ad hoc 
low-energy cutoff. Rather we encourage the computation of the injected electron spectrum 
through a more realistic model which includes both the effects of friction and diffusion 
in the evolution of the electron distribution as a whole, letting the low-energy end of the 
injected spectrum (and hence the injected power in nonthermal electrons) be determined 
naturally from the underlying physics. 


2. RELATION BETWEEN THE SOURCE-INTEGRATED SPECTRUM 
AND THE INJECTED ELECTRON SPECTRUM 


Brown ct al. 


(120031) have shown that the bremsstrahlung hard X-ray spectrum 1(e) 


(photons cm 2 s 1 keV 1 at the Earth) is related to the source-integrated electron flux 
spectrum ( nVF)(E ) (electrons cm -2 s -1 keV -1 ) by 


2 r oo n -j roo 

/(e) = Sip/ L n ^ F ( E ’ r 'l dV ff ( £ . E)dE= — 2 j' (nVF)(E)a(e,E)dE (1) 

Here n is the ambient proton density (cm -3 ) at position r, V (cm 3 ) is the source volume, 
R — 1 AU, and a(e,E) (cm 2 keV -1 ) is the bremsstrahlung cross-section, differential in 
photon energy e. For a stratified one-dimensional target, 

(nVF) (E) = A f F(E, N) dN , (2) 

Jv 

where F(E, N) is the electron flux spectrum, A is the cross-sectional area in the direction 
of electron propagation z, and N = J n(z) dz is the column density (cm -2 ). 
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For a specified form of the bremsstrahlung cross-section a(e,E), Equation (pQ) uniquely 
determines (nVF)(E) from observations of /(e); no assumptions are required regarding the 
dynamics of the emitting electrons. On the other hand, relating (nVF)(E) to the injected 
electron flux spectrum F 0 (E 0 ) (electrons cm" 2 s” 1 keV -1 ) does require assumptions on the 


dynamics of electrons in the 


target . 


the cold-target value flBrown 


2009 


Kontar et ah 


19721: 


and any change in the energy loss (or gain) 


Emslie 


1978) will change this relation (e.g., 


com pared to 


Brown ct ah 


201211 . 


2.1. Form of the injected spectrum for a prescribed source-integrated electron 

spectrum 


Galloway et al 


020 050 have found that in a warm target both friction and diffusion due 


to Coulomb collisions affect the evolution of the ensemble of accelerated particles. Therefore 
they are both important in establishing the relationship between the source-integrated 


electron spectrum (nVF)(E) a nd t 


an environment, 


Jeffrey et ah 


le injected electron spectrum F 0 (E 0 ). To describe such 


( 201 4) used the Fokker-Planck equation 


A 


dF 

dz 


= 2 Kn 



8 E 2 d/i 


+ 


+ F 0 (E) 8(z) (3) 


Here F(E, /i, *) is the local electron flux at position z along the guiding magnetic field, 
energy E, and pitch-angle cosine /i. Fq(E) 5(z) is the source of accelerated electrons which 
are injected into the target at z = 0, and K = 27re 4 A is the collision parameter, with A 
being the Coulomb logarithm. G(u) is the Chandrasekhar function, given by 


G(u) 


erf(w) — w erf (u) 
2u 2 


( 4 ) 

















































where erf(w) = {2/y/n) f exp(—t 2 ) dt is the error function. 


Averaging Equ atio n (131) ov er pitch-angle and integrating over the emitting volume 


V — f Adz (as in 


Kontar et al. 


20141 ) gives the relationship between the source-integrated 


electron spectrum ( nVF)(E ) and the (pitch-angle averaged) injected electron spectrum 
F 0 (E 0 ), viz. 


FIF\- 2K d 

F„(E„) - — ^ 



d{nVF)(E) 1 / E 
~dE + E \kf ~ 


(nVF)(E) 


■ (5) 


\ / J E=Eq 

It should be noted that substitution of a Maxwellian at temperature T: (nVF)(E) ~ 

E exp (—E/kT) results in the right side of Equation (J5J) vanishing identically, so that 
the required injected flux (source function) F 0 (E 0 ) corresponding to a Maxwellian form 
of (nVF)(E) is zero, as it should be in a dynamical model which is solely governed by 
collisional effects. It also follows that adding a Maxwellian of arbitrary size to an inferred 
(■ nVF)(E ) has no effect on the resulting F 0 (E 0 ). 


Let us now compare the constitutive relation ([5]) between F 0 (E 0 ) and ( nVF)(E ) with 
those used by previous authors. Neglecting the second-order energy diffusion term in 
Equation (0) gives 


F (jp \ — ^ 

F 0 (S„) - — Jz 


1 e{w- 1]g 





W ) <' nVF)(E ) 


( 6 ) 


E=Eq 


This is the relation ship in the finite-temperature diffusionless “warm-target” model 


proposed by 


Emslie (2003). 


As a final simplification, we may assume that the target is “cold,” i.e., E kT ((Brown 


19711 ). In this regime, G[u) —> 1/2 u 2 (Equation ((4])), resulting in the compact expression 
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Wo) = - 


K d 


previously used by 


A 


dE 


Emslie &; Smith! ( 1984 1 and lBrown fc Emsliel ( 19881 ). 


(■ nVF)(E ) 
E 


(7) 


- E = En 


Although the relationship ([7]) has been used by a number of authors to deduce the 
injected flux spectrum F 0 (E 0 ) from the source-integrated electron spectrum ( nVF )(E), we 
stress that the form (f5l) is a much more accurate representation of the t rue r elationship 


in a warm target. While the relationship (jUJ), first suggested by 1 Emsl ie! (120031) . does take 


into account the effect of finite target temperature on the energy loss rate, it neglects the 
concomitant energy diffusion responsible for thermalization in a warm target, which, as we 
shall see, is a finite temperature effect of even greater importance. 


In order to evaluate the impact of energy diffusion on the relationship between the 
injected flux spectrum F 0 (E 0 ) and the source-integrated electron spectrum (nVF)(E), 
we performed the following simple analysis. First, we assumed a hard X-ray spectrum 
consisting of combination of a Maxwellian with emission measure EM = 1 x 10 48 cm" 3 and 
temperature T = 20 MK (kT ~ 1.75 keV) and a nonthermal spectrum with power-law form 
(nVF)(E) = 10 54 (i7/10 keV) y-E 1 > E*, corresponding to a moderately large solar flare. 
We selected ( = 3 and E* = 3 keV (see Figure [1]); as we shall show, the results are insensitive 
to the value of A* as long as it is significantly less than about 10 keV. Then, using a model 
warm target with kT = 1.75 keV, we used Equation (J5J) to determine the corresponding 
form of the areally-integrated injected spectrum AFq(Eq) (electrons s _1 keV” 1 ). We then 
compared this with the forms of AF 0 (E 0 ) resulting from application of Equations (lf)l) 
(non-diffusional warm target) and ([T]) (cold target), respectively. 

Figure [T| shows the results, with the salient features summarized below. 


1. At high energies, the cold target result (solid line) is, as may be verified analytically 














10 - 



Energy [keV] 



Energy [keV] 


Fig. 1.— Top: Mean electron flux prescribed by ( nVF)(E) = 10 48 F M (E) + 
10 54 (-E/10 keV)~ 3 (electrons s” 1 keV^ 1 ), where Fm{E) is the Maxwellian flux spectrum 
normalized to unity with kT = 1.75 keV. Bottom: Forms of the injected spectrum AF 0 (E 0 ) 
recovered for a prescribed top-panel (nVF)(E) using different model assumptions: solid line 
- cold target (Equation (J7J) ), dotted-line: non-diffusional warm target (Equation (EJ), and 
dashed line: diffusional warm target (Equation (J5J) ). 
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from Equation (JTJ) , a power law with index S = ( + 2 = 5. Extending this spectrum 
down to the 3 keV lower boundary of the (nVF)(E) spectrum requires an injection 
rate N = A / 3 °^ eV F 0 (E 0 ) dE 0 ~ 10 3 ' electrons s _1 and a corresponding power 

V = A J 3 °^ cV E 0 F 0 (E 0 ) dE 0 ~ 7 x 10 28 erg s _1 . Because of this steep spectral form, 
imposition of an arbitrary low-energy cutoff E c is necessary to keep the injected power 

V = A E 0 F 0 (E 0 ) dE 0 at an acceptably small value. 


2. The result for a diffusionless warm target (Equation (El) ) do es fla 


keV and imposes a lower bound to E c at E 0 ~ kT (lEmslie 

Figure ID)- 


2003 


den off below a few 
bottom panel of 


3. Application of the diffusional warm target model (Equation [5]) causes the Maxwellian 
component in the hard X-ray spectrum to have a corresponding null injected 
spectrum F 0 (E 0 ); only the nonthermal power-law component is associated with 
a finite F 0 (E 0 ). Further, the inferred F 0 (E 0 ) has a rather abrupt cutoff at 
Eq = E C)e s ~ 9 keV, and the F 0 (E 0 ) is also noticeably flatter than the cold-target 
result up to energies as high as twice this cutoff value. Because of these features, 
this model requires only N = A F 0 (E 0 ) dE 0 ~ 2 x 10 35 electrons s~ 3 , a factor of 
~50 lower than in a cold target model above 3 keV. Similarly, the required power 
V = A E 0 F 0 (E 0 ) dE 0 ~ 4 x 10 27 erg s _1 is a factor of 20 lower than the power 
above the a priori imposed cutoff energy of 3 keV. We note that the form of F 0 (E 0 ) 
allows the integrals for the injected rate and power to be taken over the entire E$ 
range [0, oo). 


It must be stressed that the effective cutoff in the form of Fq(E 0 ) at Eq = E C)S s ~ 9 keV 
results naturally from the physics, rather than being imposed a posteriori in order to keep 


the injected power at a minimum level (e.g. maximum value of the cut-off, 


Holman et al. 
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2003 ). To see this, we assume a source-integrated mean electron fluoQ ( nVF)(E) oc E “q 


J 


and set the right hand side of Equation (j5J) to zero: 


d 

dE 



Gl A 


^d + ifT_ 1)s -c 

dE E l kT 1 


= 0 . 


E=Eq 


For E kT , G (\jE/kT\ ~ kT f‘IE and we can simplify further: 


( 8 ) 


— U-c-2 /A_( C + i)j 
dE [ \kT 1 J JJ 


= 0 , 


E=E 0 


which vanishes at an effective cutoff energy Eq = E c ^ e g, where 


( 9 ) 


Bo,.II = « + 2 )kT . (10) 

For kT = 1.75 keV and ( = 3, Equation (ITU gives E c>e s — 9 keV, in agreement with 
the results shown in Figure Q] The effective low energy cut-off E c eS could be as large as 
24 keV for a relatively soft spectrum with ( = 6 and a relatively hot (kT ~ 3 keV) flaring 
plasma; such a value of E r Bf f is comparable to the typical (upper bound) values used in the 


literature (e.g., 


Holman et al. 


2011 !). 


We further note that below the effective cutoff energy E Cie ff, the values of F 0 (E 0 ) 
inferred from Equation (J5J) are negative. This unphysical situation is forced by the 
assumption of a strict power-law form of ( nVF)(E ) down to low energies E*. In a real flare 
(with positive semidefinite Eo(Eo)), the form of ( nVF)(E) will necessarily deviate from this 
power-law form and instead transition to a Maxwellian form at low energies. This will be 


2 Note that adding a Maxwellian with temperature kT and arbitrary magnitude does not 
affect the result. So one can assume that we have a Maxwellian plus power-law mean flux 
spectrum, as is often observed in solar flares. 
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studied in the following subsection. 


2.2. Form of the source-integrated electron spectrum for a prescribed injected 

spectrum 

Considerable insight into the effects of the energy diffusion term can be obtained by 
carrying out the reverse procedure to the previous section, i.e., solving Equation (JSJ) for 
(nVF)(E) for a given injected spectrum F 0 (E 0 ), and comparing the results with numerical 
simulations based on the full Fokker-Planck equation (j3J). 

A first integral of Equation ([5]) is 


d(nVF)(E) 

dE 


1 1 
Tf~~e 


(nVF)(E) = 


A 


2 KG 


Fq(E 0 ) dE 0 


( 11 ) 


Adding an integrating factor e E//fcT allows this to be written 


d_ 

dF 


0 E/kT 

E 


(nVF) (E) ) = ( — 


0 E/kT 


EG 


AFq(Eq) dEo , 


( 12 ) 


which can in turn be straightforwardly integrated from a lower limit E m ; n (see next 
subsection) to E to give the desired expression: 


(■ nVF)(E ) 


— Ee~ E ' kT 

2 K 


e E '! kT dE' 


E'G 



A Fo(Eq) dEo . 


(13) 


It is important to notice that for E min = 0 the solution (fl3l) is infinite ; this can be 
seen by considering the limit for E' —* 0 and noting that G{u) —> u as u —* 0, so that the 
outer integrand ~ E '~ 3 ^ 2 and hence the integral diverges as E m ; n —> 0. Physically, this 
divergence arises from the fact that a finite stationary solution does not exist because we 
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have a constant source of particles, no sink of particles, and a Fokker-Planck collisional 
operator that conserves the total number of particles. Hence, in the absence of escape, the 
number of electrons in the stationary state grows indefinitely. However, as we shall now 
show, the finite extent of the warm plasma region leads to an escape term which results in 
a non-zero lower limit E m \ n and hence a finite number of electrons. 


2.2.1. Effects of escape from a finite-extent hot plasma region 


In a standard flare scenario, the electrons accelerated in the hot corona propagate 
down wards towards the d ense chromospheric layers. Recent high resolution observations 


(e.g., 


Battaglia & Komar 


20111) confirm this picture and show that the ~40 keV electrons 


form HXR footpoints (height ~1 Mm above the photosphere), below the position of the 
EUV emission (height ~(2-3) Mm) produced as a result of energy deposition by lower 


energy electrons or therma 
a number of authors (e.g., 


conduction. From a theoretical standpo int, as shown by 


Nagai fe Ernslie 


1984 : 


Mariska et ah 


19891) . the injection of 


electrons into the top layers of the chromosphere could result in a rapid rise of temperature 
of the preflare chromospheric material from T ~ 10 4 K (kT ~ 0.001 keV) through the 


region of radiative instability (ICox fe Tucker 


1969|), to a temperature] T ~ (1 — 3) x 10 7 K 


(kT ~ (1 — 3) keV) in a few seconds or slower depending on the power injected by the beam. 
However, below a certain height, the density of the chromosphere is so high that the plasma 
can effectively radiate away the injected electron power. Therefore, both observations and 
theoretical expectations show that any realistic flare volume will be composed of two main 
parts: a portion of the target (low chromosphere) that is cold {kT ~ 0.001 keV <C E ) and 


3 The flare coronal temperature is effectively limited to a few xlO' K because of the very 


strong (oc T'/ 2 ; 


Spitzer 


1962) dependence of thermal conductive losses on temperature. 
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another portion that is “warm” (E ~ kT). In the latter volume the effect of diffusional 
thermalization of the fast electrons is important. 


Therefore let us consider the warm target portion of the loop, with density n and finite 
length L and considered sufficiently “thick” to thermalize the low-energy region of the 
injected suprathermal particles, while higher-energy electrons may still propagate through 
this warm target without significant energy loss or scattering. The spatial transport of such 
electrons in the warm targ et is dominated by collisional diffusion (similarly to Equation 


(22) in 


Kontar et al. 


2014]), giving 





. . E \ dF 1 ( E 

X kT j dE + E [kT 




- F 0 (E)5(z ) , 
(14) 

where D c is the collisional diffusion coefficient. Since the collisional operator, the second 
term in Equation (ITT]) , conserves the number of electrons, the effect of injection of electrons 
into such a coronal region at a rate N (s~ x ) is to increase N, the number of electrons in 
the target. I 11 a stationary state the number of electrons in the target is balanced between 
injection and diffusive escape of thermalized electrons with energy E ~ kT out of the warm 
part of the target, at a rate 


where 


N 


d_ 

dz 




2 kT r e 
L 2 m e 


N , 


(15) 


T P = 3 


lnm e ( kT ) 3 / 2 


is the electron collision time (e.g. 


Lifshitz & Pitaevskii 

1981; 

Book 

1983) 


(16) 


19831 ). This balance 


means that both N and the source-integrated electron spectrum ( nVF)(E ) are now hnite, 
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corresponding to a finite value of E min in Equation (1T3|) . 

To estimate the value of E min , we proceed as follows. At thermal energies E ~ kT the 
approximate solution of Equation (TT3|) is 


{nvFm *w JlE Ee ~ E,tT - W-^ Ee ' mT ^ (i7) 

which has a Maxwellian form. Comparing this to the normalized expression for a Maxwellian 
with number density n (cnr 3 ) and temperature T, viz. 

<nl/F)(£) = \/Sp^ s ^ /lT ' (18) 

we see that N, the total number of thermalized electrons in the target, and N, the injection 
rate of non-thermal electrons into the target, are related by 


2 K V E, 


3y/n I kT ^ I 8 nN 


7T m e 


(kT) 3 / 2 


or 


Y _ 2 / E m[n N 


(19) 


Ek V kT r e 

Finally, comparing the diffusion result (Equation f(l5jl ) with Equation f(l9jl . we obtain an 
approximate expression for E min : 


kT L 2 m 2 

This can be written in the equivalent form 


\y/2KnL ) 


( 20 ) 


kT 


5A 


( 21 ) 
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where A = (kT) 2 /2Kn, the distance required to stop an electron of energy kT in a cold 
target. This quantity is proportional to the collisional mean free path, and we shall 
henceforth refer to it as the “mean free path.” 

The value of E m in is very sensitive to the target temperature T ( E min ~ T 9 ). ft also 
depends inversely as the fourth power of the column density £ = nL in the high temperature 
(coronal) region of the target. This shows that, unlike in the cold-target model (in which 
the relationship between F 0 (E 0 ) and ( nVF)(E ) is independent of the spatial structure of 
the source - see Equation (J7J) ), the extent L of the warm target region directly influences 
the relationship between the injected electron spectrum F 0 (E 0 ) and the source-integrated 
electron spectrum ( nVF)(E ). In particular, for sufficiently extended coronal sources, E min 
is very small and so from Equation (fT3l) . ( nVF)(E ) is greatly enhanced compared to its 
cold-target value (see Figure Ws- 

The effect of thermalization is particularly important when the hot plasma region can 
stop a large fraction of the injected electrons, i.e., when E 2 <2KnL ( E c being the lower 
cutoff energy in the injected electron spectrum; see Equation (1231) below). In the opposite 
limit, E 2 > 2 KnL, injected electrons will reach the cold chromosphere without signiheant 
thermalization. For a large injection rate N, the thermalization of injected electrons could 
be so significant such as to directly increase the density of the target. Specifically, if the 
deduced value of N satisfies 

7 £ iV > J_ [E~ 2 /5AA 2 

nAL ~ -*/¥ V kT yJT \ L ) ’ 

the plasma heating by the injected electrons needs to be taken into account. In other words, 
the minimum value of the low-energy limit E m[n is 
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Fig. 2.— The spatially integrated electron spectra ( nVF)(E ) given by the solution (ITBli 
for various values of E min given by Equation (liZTj) , for a temperature and number density of 
T = 20 MK and n = 10 11 cm -3 , respectively. The injected spectrum AF 0 (E 0 ) is a power law 
with 5 = 4 and an injection rate N = 10 36 s _1 above 10 keV. The cold thick-target model is 
represented by the light blue line, and the vertical line at 3 keV indicates the lower limit of 
the RHESSI observation range. 
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E • 7r 

-^min n 

~kT ~ 4 



where AL is the volume of hot plasma. 


( 22 ) 


2.3. Numerical Simulations 


In order to ascertain the accuracy of th e ana lytic results above, we performed 


simulations using the methodology of 


Jeffrey et al . 


(2014), which includes the effects of 


both energy and pitch angle diffusion (see Equation f[3]) ). Mimicking a real solar flare, we 
considered contiguous regions containing warm (coronal and/or flare-heated chromospheric) 
and cold (chromospheric) plasmas, respectively (see left panel in Figure [3]): 


1. Cold plasma at \z\ > L, the temperature T falls to a low T < 1 MK and the density 
rises to 1 x 10 12 cm -3 ; 

2. Hot plasma at \z\ < L , where warm target effects are important. 


Beyond \z\ = L the backgroimd represents a region with properties similar to that of 
the lower chromosphere during a flare, with a temperature T much lower than that of the 
flaring corona and a much higher number density n, that can collisionally stop electrons in 
the energy range in question over a very short distance -C 1". In this region, electrons with 
energies less than 1 keV are removed from the simulation, since at such low T and high n, 
it is unlikely that they will find themselves back in the coronal region between — L < z < L. 

The resulting stationary solution is then given by the mean electron flux spectrum 
(■ nVF(E )) = AY^nF^E) A z, where A z is the bin size and A is the cross-sectional area of 
the loop. The sum is over all simulation steps, for all z and /i (including electrons in the 
chromospheric region \z\ > L with E > 1 keV). 
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2.3.1. Simulation results 

For all simulation runs, we injected 10 4 electrons, with an injected electron spectrum 
of the form 


0 


Fq(E 0 ) — < 


CEq S 


(23) 


; Eq < E c 
; E c ft E 0 ft E max 
0 5 E 0 t> E m ax 

with spectral index 5 = 4, and low and high energy cutoffs of 10 and 50 keV, respectively. 
As a reference calculation, we used a coronal background temperature of T = 20 MK and 
a background number density n = lx 10 11 cm -3 within a region \z\ < L — 20", so that 
the actual length of the region within the bounds is 40". This is large enough to contain 
the thick targe t coronal X-ray sources of 30" or so as seen in RHESSI observations (e.g., 


Xu et al. 


20081 ). The temperature distribution T(z) and number density distibution n(z) 


over the entir e region is s’ 


method see 


Jeffrey et al 


rown in the left panel of Figure [3] For details of the simulation 


fl2014h . 
















-40 -20 0 20 40 

z [orcsec] 



Energy [keV] 



Fig. 3.— The source-integrated electron spectra (nVF)(E) for a simulation with T = 20 MK, n = lx 10 11 cm” 3 , 
6 = 4, an injected electron source size of 1" and an injection rate N — 1 s” 1 . Left-. The model temperature T(z) 
and number density n(z) distributions. Middle : No boundary conditions are imposed and hence the total integrated 


number of electrons continuously grows as the number of simulation steps increases. Hence, a stationary solution 


cannot be obtained. Right: Using the same simulation, but now applying the boundary conditions listed in Section 


12.31 a stationary solution is now produced. 


The numbers in the middle and right panel legends denote the spatially 


integrated spectrum at the given simulation “step.” 
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Fig. 4.— The spatial distribution of X-ray emission along the coronal and chromospheric 
parts of the loop for two plasma densities of 10 10 cm -3 and 10 11 cm -3 . The grey area 
delineates the chromospheric part of the loop. The parameters of the injected electrons are 
the same as in Figure |3] 

The middle panel of Figure [3] shows a simulation with no boundary conditions imposed; 
the ever-increasing growth of the number of particles in the target is evident. By contrast, 
the right panel of Figure [3] shows the simulation results with the boundary conditions above 
imposed. The electron flux distribution reaches a stationary solution long before the end of 
the simulation. Figure 0] shows the resulting spatial distribution for two different simulation 
runs where the only difference is the coronal number density in the region \z\ < L. Lower 
coronal densities lead to stronger footpoint emission. 


2.3.2. Effect of varying the low-energy cutoff 











Fig. 5.— Left panel. Two different injected electron source distributions F 0 (E 0 ): (1) a power law with 5 = 4 and 
E c = 3 keV (dashed), and (2) a power law with <5 = 4 and E c = 10 keV (solid). Both have the same total number 
of injected electrons N — lx 10 36 electrons s _1 . Right panel: the resulting source-integrated electron flux spectra 
(nVF)(E), obtained from Equations (fT3]i and (T2T|) . for the case T = 20 MK, n = lx 10 11 cm 3 , L = 20", source volume 
V = 10 27 cm 3 , corresponding to the background distribution shown in orange. Results for ( nVF)(E) are denoted by 
the black dashed curve (spectrum (1)) and the black solid curve (spectrum (2)). The cold target results are shown in 
blue (dashed and solid). The total distributions (black plus orange) are shown by the green (dashed and solid) curves. 
The lower limit of the RHESSI data (at E ~ 3 keV) is shown as a vertical line. 
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In the left panel of Figure [5] we show two different injected electron distributions 
F 0 (E 0 ): a truncated power law with 5 = 4, E c = 3 keV and a truncated power law with 
8 — 4, E c — 10 keV. In the right panel of Figure 0 the resulting mean electron flux 
spectra ( nVF)(E ) are calculated for both cases using Equations (1131) and fl2T| . The cold 
target results are also plotted for comparison, as is the thermal background distribution 
(nVF) t h(E) for a source volume of V = 10 27 cm 3 . The right panel of Figure 0 also shows 
the overall total ( nVF)(E ) = {nVF)^ eam {E) + ( nVF) t h{E ). While the solutions are 
significantly different at energies E> 10 keV, they are remarkably similar at lower energies, 
a consequence of the thermalization process since the electron behavior is mostly controlled 
by background Maxwellian properties. 


2.3.3. Comparison with previous results 


In Figure 
the results of 


El we 


compare the simulation results (using boundary conditions) with 


Ems lie (120031 1. where the case of a warm target without diffusion was 


studied analytically. The mean source electron flux spectrum is plotted for three target 
temperatures of T = 15,20 and 30 MK, using a number density of n = 1 x 10 11 cm -3 ; the 
cold target result is also included for comparison. The unrealistic case of a warm target 
without diffusion results in a large number of electrons accumulating at energies below 
that of the simulation thermal energy since such a model does not contain the physics to 
adequately describe the behavior of electrons E ~ kT . Once diffusion is added, this feature 
disappears, since the energy can now diffuse about this point leading to the formation of a 
thermal distribution at lower energies. Hence the resulting spectrum has a large thermal 
component that dominates even above E — 10 keV, making the resulting spectral index 
steeper between 10 and 30 keV. 
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Fig. 6.— Comparison of the spectra created from simulations with various physical charac¬ 
teristics (solid: warm target including diffusion; dashed : warm target, no diffusion; and grey 
dashed-dot: cold target). All curves are normalized for an injection rate of IV = 10 36 s _1 . 
The plasma parameters are T = 15 MK (orange), T = 20 MK (green) and T = 30 MK 
(blue), and n = 1 x 10 11 cm" 3 in all cases. The initial electron distribution has a spectral 
index 5 = 4 and an injected energy range from E c = 10 keV to E max = 50 keV. The vertical 
dashed line at 10 keV indicates the low energy cut-off. 
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Fig. 7.— Comparison of the analytic and simulation results for the mean source electron 
spectrum (nVF)(E). The background temperature is T = 20 MK, with a (half) loop length 
of 20". The electron source parameters are <5 = 4, iT c = 10 keV, E m ax = 50 keV and N = 10 36 
electrons s _1 . Each plot has a different background density, from n = lx 10 10 cm -3 ( top left), 
n = 5 x 10 10 cm -3 ( top right), n = lx 10 11 cm -3 ( bottom left) and n = 5 x 10 11 cm -3 ( bottom 
right). Pink: simulation, black: analytical result (TT3]i with E m[n given by Equation (l2T]h 
green (in applicable [high density] cases): best match to the simulation using a value of E min 
increased by a factor of 3 (see Equation (|25|i ). and light blue: cold target. 

























Energy [keV] 


Energy [keV] 


Energy [keV] 


Fig. 8.— As for Figure [71 but a comparison between the simulated and analytical results is shown for different 
background temperatures: T — 15 MK (left), T — 20 MK (middle) and T — 30 MK (right). In all cases the background 
density is n = 10 11 cm" 3 . All other parameters are as noted in the Figure [ 7 ] caption. Pink: simulation, black: analytical 
result (TT3l) with E mm given by Equation (T2TT) . green: best match to the simulation using a value of E min increased by a 
factor of 3 (see Equation (l25|l ). and light blue: cold target. 
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2.4. Comparison of analytical and numerical results 

We can now evaluate the accuracy of the analytic solution of Section 12.2.11 
(Equations (fT3|) and (|2T]1 ) by comparing it with the simulation results of Section [2.3.11 
For this purpose, three sets of simulation runs were performed. In Figure [7] we compare 
the analytical and numerical solutions for a case with temperature T = 20 MK, and for 
densities n = 10 10 cm -3 , 5 x 10 10 cm" 3 , 10 11 cm" 3 and 5 x 10 11 cm" 3 , respectively. In 
Figure |H] we repeat this comparison for a plasma number density of n = 10 11 cm" 3 and three 
different temperatures (T = 15 MK, T = 20 MK and T = 30 MK). In Figure [9] we compare 
the numerical results with simulation runs using different plasma parameters beyond the 
boundary at \z\ = L and different electron source parameters. In the loop region, we again 
have T = 20 MK and n = 10 11 cm" 3 , but 

1. we inject the electron source over a 10" region instead of 1", or 

2. we change the temperature beyond the spatial boundary L to T = 8 MK, or 

3. we inject an electron distribution that is initially completely isotropic, instead of 
initially beamed. 

Overall Figures [T] IHJ a n d [9] show that the analytic result matches the simulation results 
well, at least for all cases where the target density is high (n^3x 10 10 cm" 3 ). However, as 
expected, the agreement becomes poorer for the low density cases (e.g., n = lx 10 10 cm" 3 ) 
when the density of the loop region moves from a thick-target regime (F? <C 2 KnL) to a 
thin-target one ( E1 2KnL). To explore this further, in Figure flOl the predicted value 

of E min /kT is plotted against the mean free path A (Equation (12T| ): the E m in values that 
best match the simulation results are also shown. The E m \ n values required to almost 
perfectly match the simulation results (for the high density (n>3 x 10 10 cm" 3 ) cases) are 
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Fig. 9.— As for Figures [7] and EJ but for different electron source and boundary conditions: 
top left: (T = 20 MK, n — lx 10 11 cm -3 ) within the region, (T — 8 MK, n — lx 10 12 cm -3 ) 
outside the boundary, and an injected electron source size of 1", top right: (T = 20 MK, n = 
1 x 10 11 cm" 3 ) within the region, (T — 8 MK, n — lx 10 12 cm" 3 ) outside the boundary and 
an injected electron source size of 10", bottom left as for top right but with n = lx IO 10 cm" 3 
within the loop region, and bottom right: (T = 20 MK, n — lx 10 11 cm" 3 ) in the region, 
(T = 0.01 MK, n — lx 10 12 cm" 3 ) outside the boundary, an injected electron source size 
of 1", but using a pitch angle (/r) distribution that is initially isotropic instead of initially 
beamed. Pink: simulation, black: analytical result (fl3l) with E m i n given by Equation (12T]) . 
green (in applicable [high density] cases): best match to the simulation using a value of E min 
increased by a factor of 3 (see Equation (177D ). and light blue: cold target. 
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in general equal to ~ 3 times the value of E mm given by Equation f|2l]l : this is an acceptable 
adjustment given the approximations made in Section 12.2.11 (For low density cases, the 
value of E min must be adjusted by a much larger factor; see grey dot in Figure HUJ) Analytic 
results using this revised value of F min are shown, for the pertinent high-density cases, in 
Figures [7] through |9] 



X=T 2 /2Kn [arcsecs] 


Fig. 10.— The orange line shows the analytical-model value of E min /kT (Equation (T2TT) ) 
versus the mean free path A, while the blue dots show the E min values that produce a mean 
electron flux spectra that best matches the results of the numerical simulations. It is apparent 
that setting F min equal to approximately three times its analytic estimate (Equation (I2T|) ) 
produces a very good agreement between the analytic and numerical results. For higher 
values of A (lower coronal densities), E m i n begins to deviate much more from the analytic 
estimate (Equation (j2TT) : grey dot), as the hot region moves into the thin-target domain. 


In summary, the solution (fT3l) . repeated here: 
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(nVF){E) = -^Ee- E ' kT 


" E e E '/ kT dE' 


®min E' Gr 


S' jE ' 


A F 0 (E 0 ) dEo , 


(24) 


with an adjusted value of E mhl given by 


-^min 

~w 


~ 3 



(25) 


reproduces the results of a full Fokker-Planck simulation extremely well, at least for 
relatively high density (n> 3 x 10 10 cm -3 ) sources. It incorporates (see Figures 0 0 and [9]) 
the power-law form of the spectrum at high energies, the Maxwellian form at low energies, 
and the transition between these two limiting regimes. 


3. SUMMARY 

Our aim in this paper was to understand the impact of collisional energy diffusion, and 
hence of therma li za/tion in a target of finite temperature, on the deduction of the properties 
of accelerated electrons in solar flares. The results show that energy diffusion dramatically 
changes the energy evolution of electrons with energies E < 10 kT. 

Via numerical simulations of the electron dynamics that take into account the effects of 
deterministic energy loss and diffusion in both energy and pitch angle, we have determined 
the form of the source-integrated electron spectrum ( nVF)(E ) for a variety of assumed 
injected spectra F 0 (E 0 ) and target models. We have also derived an analytic expression 
(Equations (1241) and ([251) ) for ( nVF)(E) that is valid in the limit where most of the injected 
electrons are collisonally stopped in the warm target region, forming a thermally relaxed 
distribution which undergoes spatial diffusion while escaping to the cold chromospheric 
region. When E m[n becomes comparable to kT , our model approximates the standard 
thick-target results. Overall, the predictions of this model compare favorably with the 
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numerical results. 


As shown in Section 12.21 the use of the more accurate relation © between the injected 
electron flux spectrum Fq(Eq) and the observationally-inferred mean source electron 
spectrum ( nVF)(E ) results in an order of magnitude reduction of the deduced number 
of injected electrons at energies E 0 < 10 kT compared to the cold target result © and 
even compared to the non-diffusional warm target result ©. Use of a physically complete 
warm-target model leads to a lower bound on the value of the low-energy cut-off E c and 
hence an upper bound on the total injected power V = A f Eq Fc(E 0 ) clE 0 . This provides 


a new alternative to the current practice (e.g., 


Holman et al. 


2003 1 of identifying the 


maximum value of E c that is consistent with the observed hard X-ray spectrum (in a 
cold-target approximation), and hence the determination of an lower bound on the total 
injected power V = A E 0 F 0 (E 0 ) dE 0 . 


We therefore discourage the use of cold thick-target model, especially in cases of 
warm and relatively dense coronal sources. Instead, we advocate the use of the more 
physically complete target model, including the effect of electron thermalization. To derive 
the source-integrated electron spectrum ( nVF)(E ) (and so, from Equation ©, the hard 
X-ray spectrum /(e)) for a prescribed injected flux spectrum F 0 (E 0 ), use Equation (1241) 
in conjunction with Equation (l25lb The development of such fit model compatible with 
RHESSI software is currently underway. 


The authors are grateful to Gordon Holman and Brian Dennis for helping to improve 
the manuscript. EPK and NHB gratefully acknowledge the financial support by an STFC 
Grant. NLSJ was funded by STFC and a SUPA scholarship. AGE was supported by NASA 
Grant NNX14AK56G. 
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